A1 = readmatrix('NiI2_LD_T.csv')

x1 = A1(:,1)
y1 = A1(:,2)
z1 = A1(:,3)*2*1e-4
y1 = y1./z1+1.5e-3
hold on
plot(x1,y1,'Linewidth',1,'Color',[0 0 0.8])
%scatter(x1,y1,5,"filled","black")
xlim([28, 72])
xticks([30:10:70])
ylim([0,0.013])
yticks(0:0.004:0.02)
xticklabels({})
yticklabels({})
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
box on
hold off

set(gcf,'Position',[0, 0, 180, 180])